Constraints on the two-particle distribution function due to the permutational 
symmetry of the higher order distribution functions 
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We investigate how the range of parameters that specify the two-particle distribution function is 
restricted if we require that this function be obtained from the n**^ order distribution functions that 
are symmetric with respect to the permutation of any two particles. We consider the simple case 
when each variable in the distribution functions can take only two values. Results for all n values 
are given, including the limit of n — > oo. We use our results to obtain bounds on the allowed values 
of magnetization and magnetic susceptibility in an n particle Fermi fluid. 
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I. INTRODUCTION 

Two particle coordinate distribution function is one 
of the central objects of statistical mechanics Any 
such distribution must satisfy the following requirements. 
Firstly, it must be non-negative everywhere. Secondly, it 
must be normalized. Thirdly, for systems of n indistin- 
guishable particles, it is obtained by integrating out (or 
summing out) n — 2 coordinates of the n particle distribu- 
tion function that is invariant under the permutation of 
any of its particles. While the first two requirements give 
explicit constrains on the possible forms of two particle 
distribution functions, this is not the case for the third 
requirement. Clearly, it follows from the third require- 
ment that the two particle distribution must be symmet- 
ric with respect to the permutation of the two particles. 
However, the third requirement is more restrictive, i.e. 
not all symmetric non-negative and normalized two par- 
ticle distribution functions can be obtained from symmet- 
ric n particle distribution functions. Therefore, it is of in- 
terest to understand what are the explicit constraints on 
the possible forms the two particle distribution function 
(other than its symmetry) that follow from the third re- 
quirement. This knowledge can be useful in various ways. 
For example, in the case of equilibrium distributions it 
would give us an opportunity to tell which features of the 
two particle distribution function are due to the permu- 
tational symmetry of the Hamiltonian and which ones 
are due to its explicit form. In particular, as will be 
shown in a specific example below, symmetry imposed 
bounds on the values of certain physical quantities can 
be obtained. Another field of applications includes solv- 
ing various correlation function integral cquations'^'^ by 
iterational procedures . For such equations initial guess 
function that incorporate the symmetry information can 
provide faster convergence. 



Our goal in this paper is two-fold. Firstly, we would 
like to work out a simple example that shows explicitly 
how the two particle distribution function is restricted 
by the fact that it is obtained from the permutationally 
symmetric n particle distribution function. Secondly, we 
want to apply this result to a physical system and show 
that it leads to certain bounds on the allowed values of 
the physical properties for this system. 



II. OBTAINING CONSTRAINTS ON THE 
TWO-PARTICLE DISTRIBUTION FUNCTION 

Consider the n variable joint probability distribution 
function fnixi, X2, Xn) in which each variable can take 
only two values, —1 and 1. The function fnixi,X2, ■■■,Xn) 
is assumed to be symmetric with respect to the permu- 
tation of its variables, non- negative everywhere and nor- 
malized. We will use the normalization in which /„ is 
normalized to one. Such distributions can describe quite 
different physical situations. For example, as discussed 
in more detail below, after suitable rescaling /„ can be 
the joint probability distribution of the z spin compo- 
nents of n particles in a Fermi fluid. It can also describe 
outcomes of a coin-tossing experiment involving n iden- 
tical coins, or its analogues. Since we will be dealing 
with reducing higher order functions to lower order ones 
it is useful to introduce parametrizations of such joint 
probability distributions that will have simple relations 
for the functions of different orders. For our purposes a 
convenient set of parameters is given by the expansion 
coefficients of fn{xi,X2,---,Xn) in terms of products of 
functions (j)Q{xi) and <j)i{xi), that are defined as follows, 
4'oixi) = 1 and (j)i{xi) = Xi. The expansion of a normal- 
ized fn{xi,X2, ■■■,Xn) has the following forrrJ. 
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where the expansion coefficients are the moments of 
fn(xi,X2, Xn) of different orders. For f„(xi,X2, Xn) 
with a given n, the series terminates at the term in- 
volving the n^^ order moment. For the functions that 
are symmetric with respect to the permutation of any 
two of their variables, moments of the same order 
must be the same. Thus, the normalized and sym- 
metric function /„(xi, 0:2, a;„) is completely defined 
by n parameters, (xi), {X1X2), {xiX2...Xn)- The al- 
lowed range of these parameters is obtained from the 
requirement that /^(xi, 2:2, > for all coordi- 
nate values. Since there are 2" such values this leads 
to 2" inequalities. However, due to the symmetry of 
fn{xi,X2, ...jXri) Only n + 1 of these inequalities are in- 
dependent. As will be shown in the specific examples 
below, these inequalities specify a closed region of the 
allowed values in the n dimensional space of parameters 

(Xi), {X1X2), {xiX2.:Xn)- 

Let us consider the allowed range of parameters for 
the function f2{xi,X2) provided that it is normalized, 
symmetric and non- negative. The function f2ixi,X2) is 
completely specified by (xi) and {xiX2). The require- 
ment of the non-negativity leads to the following three 
inequalities 

l + 2(a;i) + (a;iX2} > 0, 
l-2(a;i) + (xiX2} > 0, 

l-(xiX2) > 0. (2) 

If these relations are treated as the equalities, then they 
define three nonparallcl straight lines in the {xi), {X1X2) 
plane. The inequalities define a region inside a triangle 
whose vertices are the intersection points of the three 
straight lines. The (xi), {X1X2) coordinates for the ver- 
tices of the triangle are (-1, 1), (0, -1), (1, 1). The region 
of allowed values is shown on Fig. ^ 

Now we would like to find out how this region is re- 
duced by the additional requirement that the function 
/2(xi, X2) is obtained by summing out variable X3 in the 
distribution function fj{xi, X2-, x^) that is non-negativity, 
normalized to one, and symmetric with respect to the 
permutation of any of its three variables. Such function 
/3(xi, X2, xa) can be completely specified by parameters 
(xi), {X1X2), and (a;iX2a;3). The requirement of the non- 
negativity for f3{xi,X2,X'i) leads to the following four 
inequalities 

1 -I- 3(a;i} 3(xi2;2) -I- (2:1X2X3) > 0, 

1 - 3(xi} -t- 3(xiX2) - (X1X2X3) > 0, 

1 -f (xi) - (X1X2) - (X1X2X3) > 0, 

1 - (xi) - (X1X2) + (X1X2X3) > 0. (3) 

These four inequalities define a tetrahedral region in the 
(xi), (X1X2}, (X1X2X3) space. The vertices of the tetra- 
hedron have the following (xi), (X1X2}, (X1X2X3} coordi- 
nates, (-1, 1, -1), ( -1/3, -1/3, 1), (1/3, -1/3, -1), (1, 1, 
!)• 

If /2(xi, X2) is obtained from /3(xi, X2, X3) then the al- 
lowed range of its parameters (xi) and (X1X2) must be the 
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FIG. 1: The regions of allowed values for the /2 parameters 
{a;i) and {x-iX2)- The region inside the big isosceles triangle 
is allowed if /2 is required to be non-negative, normalized, 
and permutationally symmetric. The areas of this triangle 
above different lines represent the allowed parameter regions 
if /2 is obtained from a permutationally symmetric /„. The 
horizontal solid line corresponds to n = 3, the gray line to 
n = 4, the dashed line to n = 5, the dotted parabola to 
n 00. 



same that these parameters have in /3(xi, X2, X3). This 
latter range is given by the projection of the tetrahedral 
parameter region for /3(xi,X2,X3) onto the (xi), (X1X2) 
plane. Since tetrahedron is a convex polytope the projec- 
tion that we seek is either a convex tetragonal or triangu- 
lar region. If we plot the actual values of the (xi), (X1X2) 
coordinates of the vertices of the tetrahedron given above 
then the tetragonal region is obtained (Fig. P). Thus, the 
requirement that /2(xi,X2) is obtained from the permu- 
tationally symmetric /3(xi,X2,X3) reduces the range of 
allowed values of (xi) and (X1X2). 

We can apply the same procedure to obtain the re- 
strictions on the /2(xi,X2) parameters if it is obtained 
from any higher order symmetric /n(xi,X2, ...,x„). The 
allowed range of n parameters for /„(xi, X2, x„) is a 
region inside a simplex in n dimensions^. The allowed 
range of values for (xi) and (X1X2) is determined by pro- 
jecting this region on the (xi), (X1X2) plane. The projec- 
tion is a convex (n -I- l)-gonal region. The (xi), (X1X2) 
coordinates of the vertices of the (rt-l- l)-gon are obtained 
from the coordinates of the simplex vertices. These lat- 
ter coordinates are obtained from n + 1 inequalities that 
specify the simplex region. Inspection of the results for 
several small n values obtained by explicit calculations 
allows one to come up with the general formulas for the 
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coordinates of the polygon vertices for arbitrary n, 
2i — n , , {2i — n)^ — n 



{X1X2) 



n(n — 1) 



(4) 



where i is an integer taking the values from to n. Ex- 
amples for n = 4 and n = 5 are given on Fig. ^ 

The obtained results allow one to draw certain conclu- 
sions about possible changes of f2{xi, X2) in an n particle 
system. Suppose that the /2 (2^1, 2:2) parameters lie in the 
region that is cut out by going to an n -I- 1 particle sys- 
tem. Then adding an extra particle to the n particle 
system will necessarily change f2{xi,X2). If, however, 
the f2{xi,X2) parameters lie in the region that is not af- 
fected by going to an particle system then f2{xi, X2) 
may or may not change depending on the details of the 
extra particle addition. 

Eqs. Q can be viewed as a parametric form of the 
function {X1X2) — {xiX2){{xi)). The explicit form of this 
function can be obtained by expressing i through (a;i) 
from the first of Eqs. Q and substituting it into the 
equation for {xiX2)- This gives 
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It can be deduced from this equation that for each n the 
polygon that defines the region of the allowed values is 
inscribed in a parabola. 

For large systems it is of interest to investigate the 
limit of n ^ 00. It is easy to check using Eqs. (@J) that 
both (xi) and {X1X2) change continuously in the limit 
n — > 00. Taking this limit in Eq. Q we obtain 



{X1X2) = {xi)' 



(6) 



The corresponding region of allowed values is shown on 
Figffl Thus, as n — > 00, the broken lines formed by n — 2 
sides of each n-gone converge to the parabola given by 
Eq. ®. 

To investigate how fast the areas A^s of the polygons 
converge to their limit at n — > cxd we use the formula for 
the area of a polygon- to obtain 
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From the statistical mechanics standpoint this means 
that if we want to approximate the n — > 00 region by 
a finite n region then we need to consider rather large 
values of n. For example, if we require a 1% accuracy 
then n > 100 must be considered. 



III. EFFECT OF CONSTRAINTS ON THE 
CORRELATION FUNCTION 

It is customary in statistical mechanics to separate 
the two particle probability distribution into the uncor- 
related and correlated parts as 
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FIG. 2: The regions of allowed values for {x\) and the cor- 
relation parameter k. The region inside the figure formed by 
the thick solid lines is allowed if /2 is required to be non- 
negative, normalized, and permutationally symmetric. The 
areas of this figure above different lines represent the allowed 
parameter regions if /2 is obtained from a permutationally 
symmetric /„. The thin solid line corresponds to n = 3, the 
gray line to n = 4, the dashed line to n = 5, the horizontal 
dotted fine to n ^ cx3. 



where 
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2;2=±1 



]2[xi,X2) 
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is the one particle distribution function and the corre- 
lation function 52(2:1, 2:2) is defined by Eq. (jSJ. Let us 
consider how the above constraints on f2{xi,X2) affect 
(72(2:1, 2:2) and /i(2;i). The functions /i(2;i) and 5(2:1, 2:2) 
can be written in terms of parameters (2:1) and (2;i2;2) as 



fiixi) 
52(2^1,2:2) 



(l + (2:i)xi), (10) 
((xi2;2) - (xi)^) 2;i2;2. (11) 



72(2:1, 2:2) = /i(2;i)/i(2;2) +52(2:1, 2:2), 



(8) 



The parameter (2:1) (that completely specifies /i(xi)) 
can be viewed as the degree of inhomogeneity since /i(2:i) 
is constant when (2:1) = 0. The parameter k = (2:12:2) — 
(xi)^ characterizes the strength of the correlations. A 
surprising result of this analysis is that the line of the 
zero correlations on the (xi), (2:12:2) plane is the same 
as the limiting parabola given by Eq. ©. To see more 
explicitly how the correlation function is affected by the 
constraints, we redraw Fig. ^ in the (2:1), k space, as 
shown on Fig. |21 

Fig. 121 shows some interesting relations between the 
correlations and inhomogeneity for finite n. For a fixed 
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(xi) the allowed range of k is always a single segment 
whereas for a fixed k the allowed range of (xi) can be 
either a single segment or, for sufficiently negative k a 
few separate segments. For this latter case it is impossi- 
ble to go from one segment to another without changing 
the correlations. Interestingly, the correlations with the 
largest negative k require fi{xi) to be homogeneous for 
even n but inhomogeneous for odd n. For n — > oo, only 
the correlations with k > are allowed. In this limit, 
the allowed correlation functions must have 172(11 1) > 0, 
52(-l, -1) > 0, and 32(1, -1) < 0, .72(1, -1) < 

As can be seen on Figs. Q] 13 the constraints from 
the higher order distribution functions do not affect the 
maximum possible range of (a:;i}. This is to be expected 
since {xi) is the only parameter characterizing /i and 
for every physical /i we can construct a permutationally 
symmetric /„ as a product of n /I's. 

IV. BOUNDS ON THE MAGNETIZATION AND 
MAGNETIC SUSCEPTIBILITY IN SPIN 1/2 
QUANTUM FLUIDS 

As an example of application of our results to a phys- 
ical system we will investigate allowed ranges of aver- 
age z component of the total spin (S^) and its vari- 
ance ((5^)^) — {S^y in a Fermi fluid^ or Bose fluid with 
(iso)spin 1/2^. Here 

^^^E^^^^E-f (12) 

i i 

and sf is the z component of spin for particle i and 
is the corresponding Pauli matrix. The quantities 
(5^) and {{S^Y) — (S^)"^ are of interest because for a 
wide range of systems in equilibrium they are propor- 
tional to the magnetization and magnetic susceptibility, 
respectively. Note, however, that our results for (5^) and 
((S'^)^) — (5'^)^ apply to arbitrary states, not necessarily 
equilibrium ones. Since erf has the eigenvalues of ±1 we 
can see that function /„ discussed above can be viewed 
as the probability distribution of the erf eigenvalues for 
n particles with variable Xi denoting eigenvalues of . 

Before proceeding with our analysis we need to be sure 
that the quantum nature of the Fermi or Bose fluid does 
not impose additional restrictions on /„. This is indeed 
the case. We will not go into details of the proof here. 
The outline of the proof is as follows. It can be shown 
that the vertices of n dimensional simplex domains for 
the allowed /„ parameters correspond to the total den- 
sity matrices composed of the wave functions that are 
eigenstates of with a given eigenvalue. (For an n par- 
ticle system there are n -I- 1 different eigenvalues, as 
many as there are vertices.) Since the density matri- 
ces for the vertices exist and since any simplex domain 
is convex, the density matrices for the points inside the 
simplex also exist. Thus, erf's can be rigorously treated 
as classical variables as far as the joint distributions of 
their eigenvalues are concerned. 
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FIG. 3; The regions of allowed values for {S^)/n and 
(((5^)^) — (5^)^) /n^. The regions bounded from above by 
the big parabolic arc and from below by series of n smaller 
arcs drawn in different styles are allowed for the n particle 
Fermi fluid. The thick solid line corresponds to n = 2, the 
thin solid line to n = 3, the gray line to n = 4, the dashed 
line to 71 = 5, the dotted straight line to n ^ 00. 

Using Eq. ((T^ and replacing af with Xi we obtain for 
the average S'^ and its variance. 

{SI - ^n{x^) 

{iSn')-{Sr = ^{nin-l)(x,X2)~n^(x,r + n) 

(13) 

Using the allowed range of parameters for {xi) and {X1X2) 
we can obtain the allowed domains for (S^) and ((5'^)^} — 
(5^)^. The results are shown on Fig. O Qualitative 
physical interpretation of these results is straightforward. 
Each point of the allowed domain with ((S*^)^) — {S^y = 
correspond to the total density matrices composed of 
eigenstates of with the same eigenvalue, i. e. (S^) 
is a multiple of h/2. For (S^) that is not a multiple 
of h/2 ((S*^)^) — (S*^)^ is necessarily greater than zero 
(because not all spins are either up or down). Generally, 
for smaller (S^) values the allowed range of ((S"^)^) — 
(S^)^ is larger. 

Consider a Fermi fluid whose Hamiltonian involves the 
coupling to the magnetic field H of the form H„i = 
~mcS^'H. For the system in the thermal equilibrium 
the magnetization per unit volume, M is given by 

M=-^{S') (14) 
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If Hm commutes with the total Hamiltonian then the 
magnetic susceptibihty, x is given by 

Thus, with suitable rescaling, Fig. |31 describes the al- 
lowed domains of M and x for any model satisfying the 
above requirements. As can be expected, no such model 
can be diamagnetic (x is always > 0). 

V. CONCLUDING REMARKS 

The generalization of our analysis to the functions 
whose variable can take more than two values is straight- 
forward. If each variable of /2 (2^1, 2:2) is allowed to take 
I values, then /2(a;i,X2) can be completely specified by 
+ 1) — 2)/2 parameters. Parameters of the higher 
order functions can always be chosen in such a way that 



they include + 1) — 2)/2 parameters of f2{xi,X2)- 
Constraints on the latter parameters are obtained by pro- 
jecting the allowed parameter range of the higher order 
function on the f2{xi,X2) parameter subspace. 

The case of continuous coordinates is more complicated 
since infinitely many parameters are needed to specify the 
distribution functions. In this case, our approach can be 
used to investigate the finite dimensional subspaces of the 
f2{xi,X2) parameter space in some regions of interest, e. 
g., in the vicinity of functions with a certain characteristic 
wavelength. 
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